ELECTRON-PHONON DRIVEN SPIN FRUSTRATION IN MULTI-BAND 
HUBBARD MODELS: MX CHAINS AND OXIDE SUPERCONDUCTORS 
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■ We discuss the consequences of both electron- phonon and electron-electron couplings in ID and 

G^ ' 

■ 2D multi-band (Peierls-Hubbard) models. After briefly discussing various analytic limits, we focus 
on (Hartree-Fock and exact) numerical studies in the intermediate regime for both couplings, where 
unusual spin-Peierls as well as long-period, frustrated ground states are found. Doping into such 
phases or near the phase boundaries can lead to further interesting phenomena such as separation 
of spin and charge, a dopant-induced phase transition of the global (parent) phase, or real-space 

^ ■ ( "bipolaronic" ) pairing. We discuss possible experimentally observable consequences of this rich phase 

O ■ diagram for halogen-bridged, transition metal, linear chain complexes (MX chains) in ID and the oxide 

superconductors in 2D. 



'in ■ INTRODUCTION 



The MX chain class of materials [1] and oxide high-temperature superconductors (HTS) are strik- 
ing similar from the perspective of strong competitions for broken-symmetry ground states — bond- 
order- waves (BOW), charge- (CDW) or spin- (SDW) density- waves or antiferromagnetism (AF), or 
both bond- and spin- distortions (spin-Peierls, SP) — and the properties of doping- and photo-induced 
local defect states (kinks, polarons, bipolarons, excitons, breathers, etc.) within the same wide variety 
of novel ground states. It is becoming increasingly apparent that to correctly model these materials, 
not only must one take into account both electron-electron (e-e) and electron-phonon (e-p) interactions, 
but also that multi-band and/or multi-orbital effects are critically important. 

Experimentally, the remarkable tunability of the MX class, along with the ease of synthesis of 
single crystals, is a major advantage allowing systematic investigation between small and large po- 
laron regimes, the influence of dimensionality, and the competition with non-adiabatic and impu- 
rity/localization effects. Intrinsic defects (self-trapped local defect states or bags) in both HTS and MX 
chain materials are thought to be polaronic in nature, and can be created via doping or photoexcitation. 
Electron-hole asymmetry, which has been proposed as a driving mechanism for superconductivity (SC), 
has been experimentally observed in MX compounds [2]. Recent evidence in PtBr suggests bipolaron 
dissociation and polaronic trapping as temperature is varied [3]. Moreover, the model described below 



manifests long-period (LP) ( "superlattice" ) phases [4], possibly observed in recent experiments on MX 
compounds, and may be related to twinning or real space pairing in the cuprate HTS. Importantly, 

SC is believed to occur at the regions of "melting" between various broken symmetry ground states 
and correlated metals. Thus we anticipate the possibility of SC, or one dimensional (ID) precursors of 
SC, in MX materials when we are near metal-insulator boundaries, where bipolarons arc extended and 
overlapping, competing with a correlated metal state. Here also the tunability of the MX class can 
be utilized, leading to our current strategy of studying PtI and Pda;Nii_2:Br to tune into weak CDW 
and SDW regimes. Recent unusual results on PtI in high magnetic fields [5] may indicate such a ID 
precursor. 

Theoretically, MX materials are clearly a ID template [6] for the same many-body methods and 
parameter determinations one must employ in 2D cuprate and 3D bismuthate HTS models when non- 
bonding orbitals are neglected. Particularly all are hybridized, multi-band materials. The appearance 
of superlattice CDW ordering in MX materials [4] may be similar to observations in many HTS ma- 
terials. In other words including a microscopic driving force for finite scale "twinning" textures and 
integrating electronic degrees of freedom yield an effective anharmonic lattice dynamics. Indeed, in the 
model described below the superlattices arise from ordering of "bipolaronic" defects with respect to 
the period-4 CDW, corresponding to a "melting" of the broken symmetry state (c/. Bai_a;Pba;Bi03). 
Polaron pairing into bipolarons is a common ingredient of many current e-e and e-p theories of HTS 
materials. Further, both MX chains and (with M=Cu, X=0) Cu02 planes (and chains), can be 
described by essentially the same multi-band, tight-binding extended Peierls-Hubbard Hamiltonian 
(PHH), described briefiy below, although a 3/4-filled, 2-band (2B) ID model for the MX compounds 
or a CuO chain, and a 5/6-filled, 3-band 2D model for the CuO planes. This same model can be 
considered a 3/4-filled analog of the the organic conductor polyacetylene, model charge-transfer salts 
such as TTF-TCNQ, or be used to investigate neutral-ionic transitions [7]. Mathematically, one can 
consider the two orbitals to be on the same site, and thus this Hamiltonian is also related to the Kondo 
Hamiltonian used to describe heavy fermion materials. 

MULTI-BAND MODEL 

Our model ID, 2B, PHH focusing on the and p^ orbitals and including only the nearest 
neighbor interactions is [8,9] 



where (e, (3, U)i = {e, (3, U)m,x and ex=2eo. For a more detailed description of the parameters, the 
methods of solution, and general properties of this Hamiltonian, see Ref.s [8] and [9]. The results here 
were obtained by exact diagonalization (ED) of the many-body PHH on small systems. 

Tuning eo, to, Um, and Ux is essentially a ID analog of the theoretical discussions in 2D which 
focus on the p-d hybridization in HTS materials. There, as here, these parameters determine the 
stoichiometric ground state broken symmetry order, and the nature of electron or hole doping into those 
ground states. Indeed, the similarity is even stronger because the nominal band filling is essentially the 
same in both MX and HTS materials - 3/4-filling of 2 bands and 5/6-filling of 3 bands, respectively. 
Furthermore, both strongly AF and CDW stoichiometric compounds exist for both HTS and MX. We 
stress the MX class has the advantage of essentially continuous tunability between these extremes. 
Examining trends in related materials often yields insights that a single material focus might easily 
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Fig. 1. The to=0 phase diagram for (3x / Pm=Ux /Um=^- The dashed hne is the slice used for Fig. 4 
and + the point in Fig. 3. Note changing P translates into moving along a line through the origin. 

Table 1. The period-4 phases, occupancies, distortions, energies, and effective J scaling in the to^^O 
limit. Here Um,x=UM,x/eo, bm,x={PM,xT /{Keo), and (^(n)=dx(cos ^-sin ^)-dM(cos ^+sin ^) 
defines X(M) sublattice distortion order parameters dx{dM) with the first bond (X-M) being short, 
(constant or (—1)" terms just renormalize to and cq). For non-period-4 phases see Ref. [4]. For toT^O, 
the unpaired spins are antiferromagneticly correlated. 
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miss. Finally, we note that competition between e-p and e-e coupling is especially important for for 
doping states and can induce strong local LD even in a strongly magnetic background [9,10]. 

ANALYTIC LIMITING CASES 

In the LD CDW case typical of many MX chains, a first approximation is to approach the material 
as decoupled X=M=X trimers and M monomers. Indeed, PtCl chain absorptions are very close to the 
monomers and trimers in solution. Conversely, when e-e interactions are dominant and the splitting 

between the bands is large (IB limit), one expects a SP phase in ID [11], with an effective AF coupling 
between M sites given by Jmm = i\iM/^M, where Imm is the M-M hopping in the effective IB model. 
We discuss here the to=0 limit, which captures the essential features of the material in both these 
limits. Another analytically tractable limit (f/=0) was discussed in detail in Ref. [8]. 

The distortion and energies of the period-4 (P4) phases for to=a=0, listed in Table 1, are: two 
undistorted configurations with unpaired electrons on the M (MAF) or X (XAF) sites, one phase 
with both LD and unpaired spins (MIX) and 2 diamagnetic phases with the M (BOW) or X (CDW) 
sublattice distorted and a charge-density-wave on the opposite sublattice. The phase diagram for 
parameters we expect to be relevant to MI or NiX materials, is shown in Fig. 1. The phase diagram 
is more complex for Ux,Px^Um,/3m, where the hybridization-driven competition is most effective. 
The BOW (XAF) phase is only found for \Px/Pm\ {Ux/Um) larger than the case shown (though at 
io/co ~ 2, or smaller with a>0, a BOW phase is found numerically near the MAF phase boundary for 
these parameters). 



SPIN-PEIERLS AND FRUSTRATED PHASES 

When to=0, all spin excitations are isoenergetic. For to^O, AF correlations develop and one can 
treat the problem in terms of an effective spin model [12]. Without coupling between the two bands, 
the lower, X-like (for eo>0) band is full and non-magnetic, while the upper (M-like) band is 1/2-full 
with one electron per M site (when e-e repulsion is dominant). The on-site e-p coupling ft leads to a 
splitting of the upper band which competes with an AF ordering of the spins caused by the effective AF 
coupling {J mm) (as in the effective IB case). When hybridization between the two bands is allowed, 
the lower band is not completely full, and now there is an effective AF coupling between neighboring 
halide sites [Jxx, dominant when Ux is dominant) as well as M and X sites (Jmx), not present in IB 
models. In fact, when the splitting due to the e-p coupling /3 is on the order of U and eo, the AF state 
with neighboring M-X pairs singly occupied can become the ground state, as shown in Fig. 1. Note this 
implies, in contrast to the IB case, that the combination of e-e and e-p coupling in the 2B model drives, 
in addition to the non-magnetic CDW and BOW phases, three (competing) SP phases: one on the 
X-sublattice (XAF), one on the M-sublattice (MAF), and one involving MX pairs (MIX). Since Jmm, 
Jxx, and Jmx arc all AF couplings (J > 0), they obviously cannot all be simultaneously satisfied 
and the system is frustrated, as shown in Fig. 2. It is straightforward to derive from fourth-order 
perturbation theory in the large U limit an effective t, Jmx, Jmm, Jxx model, though the expressions 
for the J's are cumbersome and phase dependent [12]. The resultant to dependencies of the J's are 
listed in Table 1. When to is small and one is in a regime where only one of the J's is important, one 
can numerically check this estimate by comparing the energies of the singlet and triplet ground states 
(at fixed LD). Fig. 3 shows the to dependence is correctly predicted. Note the CDW phase has an 
entirely e-p driven AF component. 

In Fig. 4 we show the total energy and average lattice distortion for parameters similar to those 
used in the 2D model for CuO planes [10] (in hole notation, eo=to/2— f7M/4, (3=0, Ux/to=3, V/to=l, 
a'^ / Kto=2, and the M-M distance was held constant) with analogous results. For small Um, the ground 
state is a CDW, as expected [8,9]. As Um increases, the X no longer symmetrically distort (BOW in 
notation of Ref. [10]), then a LP phase develops, and finally the SDW phase expected at large Um sets 
in. In contrast to the 2D case, the SP phase is not seen, though a metastable SP regime is found. As is 
clearly seen from the figure, the "melting" of the CDW phase takes place through a LP intermediary 
[13]. SC is believed to occur at the regions of "melting" between various broken symmetry ground 
states and correlated metals. While SC is not expected in the ID materials, unusual results on PtI in 
high magnetic fields [5] may be a ID precursor. 

Upon doping, the phase diagram becomes even richer. As has also been seen in the 2D version of 
the PHH [10], doping into, e.g., the SDW phase can lead to polaronic defects with local CDW character, 
and vice-versa. An example of this for our ID model was reported in Ref. [9] for Ni parameters. Similar 
results have recently been reported in a 2D model for the HTS [10]. Further, near the phase boundary, 
we have also found regions where doping into the SDW ground state can cause a CDW defect in a 
CDW background (and vice-versa) - i.e, the dopant has altered not only it's local environment, but 
the background phase as well. While clearly finite-size effects will be important for this energy balance, 
one can easily find parameter regions where the "finite size" is as large as typical correlation lengths in 
the real materials. Additionally, in the "zero-hopping" limit where the many-body PHH is analytically 
tractable, one can show that such phase reversal upon doping can occur even in the infinite-system 
limit. Indeed, the narrow CDW phase seen near l/8th doping in the La based HTS compounds may 
be related to this phenomenon. 

The crossover from P4 CDW to P4 MAF may also be accompanied by LP superlatticc phases 
[14]. Such superlattice phases have recently been found in a 2D, 3B model of HTS Cu02 layers. 



Fig. 2 (left). Schematic energy level diagram in the strongly correlated limit showing frustration of the 
effective antiferromagnetic couplings caused by MX hybridization. 

Fig. 3 (right). The to dependence of the difference in singlet and triplet energies for the CDW, MIX, 
and MAF phases at the + in Fig. 1. For small to, this corresponds to excitations of the eflfective 
spin-Hamiltonian H = AJSiSi+j with energy AJmm in the MAF phase and 4Jmx in the MIX phase, 
whereas in the CDW phase this reflects the Peierls gap. 



Fig. 4. (a) Total energy and (b) lattice distortion amplitude as a function of Um for CDW (— ), BOW 
(- -), SP (• ■ •) SDW ( — ), and LP ( ) phases. Parameters consistent with CuO were used. 

including intersite c-p coupling. We anticipate the same behavior in our MX model when intersite 
e-p coupling is included. Such superlattices may be viewed as an ordered array of discommensuration 
defects with respect to a nearly stable commensurate order (e.g., P4). In view of the effective J's 
discussed above, it may be natural to model such states in terms of ANNNI-like models, where nearest 
and longer range couplings compete, leading to frustration and associated complexity phenomena - 
multitimescale relaxation, hysteresis, metastability, etc. In the context of the MX class, it will be 
particularly interesting to investigate materials in, or near, this crossover regime - e.g., PtI - and to 
further control the crossover with pressure, magnetic field [5], doping, impurities, etc. Indeed doping 
into this complex regime should be highly sensitive to the softness and competitions of the phases: 
This may well be an excellent regime to study pairing tendencies and metallization. 

CONCLUSIONS 

Wc stress that the ID, 2B, PHH, while simple to write down, is representative of a very large 
variety of low-D electronic materials. This variety is mirrored in the model's richness, especially in 
terms of doping near phase boundaries where novel pairing mechanisms are found (inclusion of e-p in 
a 2D 3B model leads to coexistence of CDW and SDW phases - "charge bags" - as e-p interactions 



are important only in the neighborhood of defects [10]), or a dopant-induced transition [15] of the 
global phase may exist, besides the usual plethora of doping and photoinduced non-linear excitations. 

The MX class of compounds arc uniquely important as a testing ground for many-body modeling and 
materials design strategy in strongly correlated, low-D materials (in particular the oxide HTS). Apart 
from pure and mixed-halide MX materials wc arc also beginning to explore mixed-metal (Ma;M{_^X) 
and bimetallic (MMX) systems, as well as effects of magnetic fields (especially on the weak CDW/SDW 
ground state materials). We feel that experimental investigations of the pressure dependence and high 
(magnetic) field behavior of pure and doped materials in the LD/AF cross-over regime, such as PtI, 
NiBr, or their mixed-metal or -halide analogs, will continue to yield interesting insights into the nature 
of multiband effects and the competition between e-e and e-p interactions. 
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